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Abstract 

Recently, the concept of cumulative residual entropy (CRE) has been stud¬ 
ied by many researchers in higher dimensions. In this article, we extend the 
definition of (dynamic) cumulative past entropy (DOPE), a dual measure of 
(dynamic) CRE, to bivariate setup and obtain some of its properties including 
bounds. We also look into the problem of extending DOPE for conditionally 
specified models. Several properties, including monotonicity, and bounds of 
DOPE are obtained for conditional distributions. It is shown that the pro¬ 
posed measure uniquely determines the distribution function. Moreover, we 
also propose a stochastic order based on this measure. 
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1 Introduction 


In recent years, cumulative residual entropy and cumulative past entropy have been 
considered as a new measure of information that parallel Shannon (1948) entropy. 
Let X be an absolutely continuous nonnegative random variable with survival func¬ 


tion F(-) = 1 — F(-) and probability density function (pdf) /(■). Then Shannon’s 
differential entropy is dehned as 



( 1 . 1 ) 


where ‘log’ means natural logarithm and, by convention, OlogO = 0. It measures 
the expected uncertainty contained in / about the predictability of an outcome of 
X. In spite of enormous success of Shannon entropy, the differential entropy fll.ip 
presents various deficiencies when it is used as a continuous counterpart of the classi¬ 


cal Shannon entropy for discrete random variable. Several attempts have been made 


in order to dehne possible alternative information measures. Recently, Rao et al. 
(2004) identihed some limitations of the use of fll.ip in measuring randomness of cer¬ 
tain systems and introduced an alternative measure of uncertainty called cumulative 
residual entropy (CRE) dehned as 



( 1 . 2 ) 


which relates to uncertainty on the future lifetime of a system. This measure is 
based on survival function and is particularly suitable to describe the information 
in problems related to ageing properties in the reliability theory. Motivated by the 


salient features of fll.2l) . Di Crescenzo and Longobardi (2009) proposed a dual concept 
of CRE called cumulative past entropy (CPE) dehned as 



(1.3) 


which measures information concerning past lifetime. 

Length of time during a study period has been considered as a prime variable 
of interest in many areas such as reliability, survival analysis, economics, business, 
etc. In particular, consider an item under study, then the information about the 
residual (past) lifetime is an important task in many applications. In such cases the 
information measures are functions of time, and thus they are dynamic in nature. 
Asadi and Zohrevand (2007) further studied the function obtained from fll.2p in the 
residual setup called dynamic CRE (DCRE), given by 



(1.4) 
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For a discussion on the properties and generalization of (dynamic) CRE one may 
refer to Rao (2005), Navarro et ah (2010), Abbasnejad et ah (2010), Kumar and 
Taneja (2011), Navarro et al. (2011), Sunoj and Linn (2012), Khorashadizadeh et al. 
(2013), Psarrakos and Navarro (2013), Navarro et al. (2014), Chamany and Baratpour 
(2014), among others. In analogy with fll.4p . Di Crescenzo and Longobardi (2009) 
also studied CPE for past lifetime called dynamic CPE (DCPE), defined as 



(1.5) 


For more properties, applications and recent developments of (II.5h . one may refer to 
Abbasnejad (2011) and Di Crescenzo and Longobardi (2013). 

In studying the reliability aspects of multi-component system with each component 
having a lifetime depending on the lifetimes of the other components, multivariate life 
distributions are employed. Reliability characteristics in the univariate case can be 
extended to the corresponding multivariate version. In a recent work, Rajesh et ah 
(2014) have considered extension of DCRE to bivariate setup and study its properties. 
Some bivariate distributions are also characterized there. Several generalizations to 
the concept of bivariate DCRE can be found in Sunoj and Linu (2012) and Rajesh 
et al. (2014). In various contexts the uncertainty is not necessarily related to the 
future but may refer to the past. Even though a lot of interest has been evoked on the 
bivariate extension of CRE, no works for CPE, to the best of our knowledge, till now, 
seem to have been done in higher dimension. It is to be noted that the concepts in past 
time are more appropriate than those truncated from below when the observations 
are predominantly from left tail. This shows the relevance and usefulness of studying 
CPE when uncertainty is related to the past. 

In the present paper we consider (dynamic) CPE for bivariate setup, and study its 
various properties useful in reliability modeling. The rest of the paper is arranged as 
follows. Section 2 includes the dehnition and basic properties of bivariate CPE. Some 
bounds for bivariate CPE are also obtained. In Section 3 we look into the behavior 
of dynamic CPE for conditional distributions. Several properties of the measures are 
studied along with monotonicity and some characterization theorems arising out of 
it. A stochastic order is proposed and studied which is based on the measures. 

2 Definition and properties of bivariate CPE 

The topic of measuring the information content for bivariate (multivariate) distribu¬ 
tions when their supports are truncated progressively are considered in recent past. 
A significant results in this area have been provided in Ebrahimi et al. (2007). In this 
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section we look into the problem of extending CPE defined in fll.Sp to the bivariate 
setup. Let X = {Xi,X 2 ) be a random vector of nonnegative absolutely continuous 
random variables with joint pdf f{xi,X 2 )- We may think of Xj, z = 1,2, as the 
lifetimes of the members of a group or components of a system. Then the Shannon 
differential entropy of (Xi,X 2 ) is dehned as 

poo poo 

H{Xi,X 2) = - dxi f{xi,X2)\ogf{xi,X2)dx2. ( 2 . 6 ) 

Jo Jo 

One of the main problems encountered while extending a univariate concept to higher 
dimensions is that it cannot be done in a unique way. A natural extension of CPE 
(ll.3p to the bivariate setup can be obtained from (12.bh by replacing f{xi,X 2 ) by 
F{xi,X 2 ) as given in the following dehnition. Since the past lifetime has always 
a hnite support we restrict our attention to random variables with hnite supports. 
Therefore, we assume that the support of (Xi,X 2 ) is included in (0,6i) x ( 0 , 62 ) for 
some nonnegative real values &i, & 2 - 

Definition 2.1 Let X = (Xi,X 2 ) be a nonnegative bivariate random vector admit¬ 
ting an absolutely continuous distribution function with joint pdf f{xi,X 2 ), distri¬ 
bution function F{xi,X 2 ), marginal distribution functions Fi{xi) and marginal pdfs 
fxii^i), i = 1,2. We define the bivariate CPE as 

pbi i>b2 

e{Xi,X2) = - / F{xi,X2)hgF{xi,X2)dx2dxi (2.7) 

Jo Jo 

provided the integral on the right hand side is finite. 


The following example clarihes the effectiveness of the proposed measure. 


Example 2.1 Let Xi and X 2 are random lifetimes of two components with joint pdf 


f{xi,X2) 


2, 0 < a:i < 1, 0 < X 2 < Xi 

0 , otherwise. 


Then, H{Xi,X 2 ) = - log2, H{Xi) = H{X 2 ) = 1/2 - log2 and e{Xi,X 2 ) = (1 - 
log2)/4, 7(Xi) = 2/9, s{X 2 ) = y ~ Here H{Xi,X 2 ) along with H{Xi) 

and H{X 2 ) are negative which do not make a sense whereas their CPE measures are 
positive. Most importantly, H{Xi) and H{X 2 ) are identical bute{Xi) ande{X 2 ) are 
different. □ 


Indeed, since log a; ^ a; — 1 for all a; > 0, we have 


£(Xi,X 2 ) ^ 


fbi pb2 


F{xi, X2) (1 - F{xi, X2)) dx 2 dxi, 


'0 Jo 
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where 1 — F{xi,X 2 ) = P(Xi > xi or X 2 > X 2 ), which gives the probability that at 
least one of the components will snrvive beyond the time {xi,X 2 )- In connection with 
Example Em the above right hand side expression is evaluated as 0.028 which verihes 
the proposed lower bound of bivariate CPE. 

We also recall from fl2.7p that if Xi and X 2 are independent then 


£(Wi,X 2 ) 


(*h2 


F2{x2)dX2 


^(W) + 


r-bi 


Fi{xi)dxi 


KX 2 ). 


( 2 . 8 ) 


uo J L~'o J 

The following additive property of bivariate CPE is due to Di Crescenzo and Longo- 
bardi (2009). 


Proposition 2.1 Let X = (Xi,X2) he a nonnegative bivariate random vector where 
Xi,X2 are independent random variables with supports [0, 61] and [0,62], respectively. 
Then 

£(Xi, X2) = (62 ~ T2) ^(-^1) + (^1 ~ Ti) ^(^2)5 
where fii = E{Xi) and bi are finite. In particular, if Xi and X2 both have support 
[0, b] and expectation p., then we have 


£(Xl, X 2 ) = {h-p) [£(Xi) + £(X2)] 


which shows that bivariate CPE also have an appealing property in analogy with Shan¬ 
non’s differential entropy for two-dimensional random variable e.g., H{Xi,X2) = 
H{Xi) + H{X2), if Xi and X2 are independent. 

Let us analyze the effect of linear transformations on the bivariate CPE. The proof 
is immediate from fl2.7l) . 

Proposition 2.2 Let Y = (Yi,Y2) be a nonnegative bivariate random vector where 
Yi = CiXi + di with Cj > 0 and di ^ 0 . Then 


£(W,E2)=CiC25(W,X2), 

which shows that bivariate CPE is a shift-independent measure. 

Now we show that bivariate CPE is not invariant under non-singular transformations. 


Proposition 2.3 Let Y = (1^,12) be a nonnegative bivariate random vector. If 
Yi = (ffiXi), i = 1,2 are one-to-one transformations with (pfixi) are differentiable 
functions, then 

rbi j^h2 

e{Yi,Y2) = - / F{xi,X2)\ogF{xi,X2)\J\dx2dxi, 

Jo Jo 

where .1 = -^(fi{xi)-^ip2{x2) is the Jacobian of the transformation. 
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A sharper lower bound for bivariate CPE is given in the following theorem. 


Theorem 2.1 For a nonnegative bivariate random vector X = (Xi,X 2 ) 

£(Xi,X 2 ) ^ max , 

where 


Cl = exp 


'•62 


Exi (log j F{Xi,X 2 )\\ogF{Xi,X 2 )\dx 2 j 


and C 2 = exp 


r-bl 


Ex 2 log / F{xi,X 2 )\\ogF{xi,X 2 )\dxi 


Proof: Using log-sum inequality we get 


C /xi(a:i)log^ 


/xi Cl) 


^ /o' fxAxi)dxi log 


F{xi,X 2)\ \ogF{x\,X2)\dX2\ 
fo^ /xi ixi)dxi _ 


dxi 


fo^ fo^ logF(xi,X2)ldx2dxi 

which on simplihcation reduces to 


= -log£(Xi,X2), 


s(Xi,X2)^Cie^(^^\ 


where Ci = exp 
one can obtain t 


Exi (log E{Xi, X 2 )I logE(Xi, X 2 )\dx 2 ) • Proceeding analogously 


he required result. 


□ 


If X = (Xi, X 2 ) represents the lifetimes of two components in a system where both 
the components are found failed at times ti and ^ 2 , respectively, then, the measure of 
uncertainty associated with the past lifetimes of the system, called bivariate dynamic 
CPE, is given by 


^ 2 ) 




E{xi,X2) 

E{ti,t2) 


log 


E{xi,X 2) 
E{ti, ^ 2 ) 


dx2dxi, 


(2.9) 


which can be thought of as two-dimensional extension of dynamic CPE. If Xi and 
X 2 are independent, then 


£x(ti,t2) = m2(t2)£xi{ti) -h mi(ti)ex2{t2) 

where miiti) = E {ti — Xj|Xj ^ U), are the marginal expected inactivity time (EIT) 
of the components Xj, i = 1,2. The bivariate dynamic CPE is also not invariant 
under non-singular transformations. 
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Proposition 2.4 Let Y = (Yi,F2) be a nonnegative bivariate random vector. If 
Yi = i = 1,2 are one-to-one transformations with (pi{xi) are differentiable 

functions, then 


£Y{^l{ti),(f2{t2)) = - 


In particular, if we choose g:)i{Xi) 
£Y{‘Pl{ti),(p2{t2)) = CiC2£x(ti,t2). 


ti ft2 


F{xi,X2) F(a;i,a;2) 


CiXi + di with Ci > 0 and di ^ 0 for i = 


1 , 2 then 


When we consider bivariate measures, it is necessary that the measurement on the 
basis of one component is not affected by the missing or unreliable data on the other 
component and hence it is necessary to consider component-wise CPE subject to the 
condition that both the components are found failed at some specihed times. Such a 
measure will be more reliable as the unreliable data is omitted. With this motivation, 
we now look into the behavior of dynamic CPE for conditional distributions. 


3 Conditional dynamic CPE 

Specihcation of the joint distribution through its component densities, namely marginals 
and conditionals has been a problem dealt with by many researchers in the past. It is 
well known that in general, the marginal densities cannot determine the joint density 
uniquely unless the variables are independent. Apart from the marginal distribution 
of Xi and the conditional distribution of Xj given Xi = ti, i = 1,2, i j, from which 
the joint distribution can always be found, the other quantities that are of relevance 
to the problem are (a) the two conditional distributions of Xi given that Xj < tj, 
i,j = 1,2, i ^ j, (b) marginal and conditional distributions of the same component 
viz. Xi and the Xi given X 2 = ^2 or X 2 and the X 2 given Xi = ti. Characteriza¬ 
tion of the bivariate density given the forms of the marginal density of Xi {X2) and 
the conditional density of Xi given X 2 = t 2 {X 2 given Xi = ti) for certain classes 
of distributions, have been considered by Seshadri and Path (1964), Nair and Nair 
(1988), Hitha and Nair (1991), Arnold et al. (2001) and Navarro and Sarabia (2013). 
Accordingly in the following Subsections 3.1 and 3.2, we consider conditional dynamic 
CPE of Xi given Xj < tj and Xi given Xj = tj, i,j = 1,2, i 7 ^ j, respectively and 
study some characteristics relationships in the context of reliability modeling. 

3.1 Conditional dynamic CPE for Xj given Xj < tj 

Here we look into the behavior of dynamic CPE for conditional distributions. Con¬ 
sider the random variables Yi = (Xj|Xi < ti,X 2 < t 2 ), i = 1,2, which correspond 
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to the conditional distributions of Xi subject to the condition that failure of the 
hrst component had occurred in (0,ti) and the second has failed before ^ 2 - The 
distribution functions of Yi are given by P ^Yi ^ yi^ = ^ < Vi < and 

P (^2 ^ 1 / 2 ^ = ’ 0 < 1/2 < h- Then the dynamic CPE for Yj, called condi¬ 

tional dynamic CPE (CDCPE), takes the form 


^ 1 ) ^ 2 ) 



-E(fi, ^ 2 ) 


log 


F{ti, t2) 


dxi 


(3.10) 


and el{X;ti,t 2 ) 



F(ti,X2) 

F{ti, t2) 


log 


F{ti,X2) 
F{ti, t2) 


dX2. 


(3.11) 


X = (Xi, X 2 ) represents a bivariate random vector, recalling (II.3p . then el{X;ti,t 2 ) 
identifies with the CPE of {Xi\Xi < ti,X 2 < ^ 2 ), with a similar interpretation 
for e 2 {X;ti,t 2 ). In particular, if Xi and X 2 are independent, then e*{X-,ti,t 2 ) = 
exi{ti), i = 1,2. In the sequel we give the dehnitions of bivariate reversed hazard 
rate and bivariate EIT functions. For more details one may refer to Roy (2002) and 
Nair and Asha (2008). 


Definition 3.1 For a random vector X = {Xi,X 2 ) with distribution functions F {ti, 12 ) 

(i) the bivariate reversed hazard rate is defined as a vector, (fulfil, t 2 ) = t 2 ), ^^(ti, 12 )) 

where (j)fi{ti,t 2 ) = ^ logF(ti, 12 ), f = 1, 2 are the components of bivariate reversed 
hazard rate; 

(a) the bivariate EIT is defined by the vector {ti,t 2 ) = {ti,t 2 ),m 2 iti,t 2 )) 

where mf {ti,t 2 ) = E (fii — Xi\Xi < ti, X 2 < t 2 ), i = l,2. Fori = l, 

1 

mf{ti,t2) = F{xi,t2)dxi, 

X [ti, t2) Jo 

which measures the expected waiting time of the first component conditioned on the 
fact that both the components were failed before times ti and t 2 , respectively. 


Note that, fl3.10l) can alternatively be written as 

/■‘i F{xi,t2) 


el{X-,ti,t2) = mf{ti,t2) logF(ti,f2) 


Jo F{ti,t2) 


\ogF{xi,t2)dxi, (3.12) 


where the above right hand side integral can be found to have a nice probabilistic 
meaning as follows. Let us set, for 0 ^ a < 6, 

rb 


T[^\a,b-,t2) ■.= - / \ogF{xi,t2)dxi. 














Its partial derivative is closely related to the distribution function of (Xi, X2). Indeed, 
from above we have = — logF(fi,^2)- Then, 


F{xi,t2) 

F{ti,t2) 


\ogF{xi,t2)dxi = 


F{ti,t2) Jo 

1 


f*Xi 


F{tl,t 2 ) Jo 
d2) 


f{u,t 2 )duj \ogF{xi,t 2 )dxi 
f{u,t 2 ) ( [ \ogF{xi,t 2 )dxi^ du 


= E^[^\X^M-M)\Xi<F,X2<t2 

So, from 03.121) it can be written that 

^1) ^2) = (^1; ^2) log F{tl,t 2 ) + E ^(Xi, ti; ^2)!-^! < tl, X 2 < t 2 


Similarly, (Id.llh can also be written as 


el{X;ti, ^ 2 ) 


mf(ti,t2)logF(ti,f2) 



T'(ti,3;2) 

F{ti,t2) 


logF{ti,X2)dX2 


(fi,t2) logF(ti,t2) + E 


Tf (X2,f2;tl)|Xi <ti,X 2 <t 2 
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where T2^^(a, 6;ti) = —logF(fi, a;2)(ix2. Differentiating 03.10p and 03.111) with re¬ 
spect to tl and ^2, respectively we get in general 

■^e*{X]ti,t2) = (j)f{ti,t2) [mf{ti,t2)-e*{X]ti,t2)] , i = l,2. (3.13) 

Now we have the following theorem. 


Theorem 3.1 Forti,t 2 > 0, e*{X;ti,t 2 ) is increasing in U, if and only if 

£*(X;ti,t2) ^ (fi,t2), f = l,2. 


Example 3.1 Let X follow the distribution as given in Example \2.1[ Then, for 
i = 1,2, mf {ti,t 2 ) = ti/2 and e*{X;ti,t 2 ) = ti/4:. Here e*{X-,ti,t 2 ) is increasing in 
ti and satisfy the above ineguality. 


Remark 3.1 Di Crescenzo and Longobardi (2009) pointed out thatex{t), defined in 
/ li.,51) . cannot be decreasing in t for any random variable X with support (0, b) with b 
finite or infinite. So, for i = I, if t 2 can be fixed at 62, then el{X;ti,t 2 ) = Sxiiti) 
cannot be decreasing in ti. Thus for all tj > 0, it can be written that £*(X; ti, ^2), i = 
1,2, i ^ j, cannot be decreasing in ti. 


Although e*{X-,ti,t 2 ) is not decreasing in ti for all tj > 0, i,j = 1,2, i ^ j, the 
following example shows that e*{X;ti,t 2 ) can be increasing in ti for all tj > 0. 
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Example 3.2 If F{xi,X 2 ) denotes joint distribution function of the random vector 
X = (Xi,X 2 ), which follows bivariate extreme value distribution of type B, then for 
all xi, X 2 ^ 0 , 


F{xi,X 2) = exp 




_j_ ^—mx2 \ 


m > 1. 


For properties of this distribution one may refer to Kotz et al. (2000). From Figure 
[H it is clear that e*{X;ti, ^ 2 ) for this distribution taking m = 2, is increasing in ti for 
all tj > 0, i,j = 1, 2, i ^ j ■ It is to be mentioned here that while plotting curves, the 
substitutions ti = —\nx and t 2 = — In y have been used so that e\{X-, — Inx, — Iny) = 
oi(x, y) and e^iX; — Inx, — In y) = 02 {x,y), say. 




Graph of ai(x, y) Graph of 02 (x, y) 

Figure 1: Graphical representations of ( 01 , 02 ) (Example I3.2p 

In the following theorem we obtain a functional relationship between the vector 
GDGPE and vector valued EIT. This relationship is useful in the sense that in many 
statistical models one may have information about bivariate EIT. 

Theorem 3.2 Let X = (Xi,X 2 ) be an absolutely continuous nonnegative bivariate 
random vector with finite e*{X;ti,t 2 ) and bivariate EIT components mf{ti,t 2 ), i = 
1 ,2. Then for all ti,t 2 ^ 0, 

mf {xi,t2)fi{xi;ti,t2)dxi 
f*t2 

and e 2 {X-,ti,t 2 ) = I (ti,X 2 )/ 2 (a: 2 Gi,G)dx 2 , 

Jo 

where fi{xi]ti,t 2 ) is the density function of {XfXi < ti,X 2 < ^ 2 ), i = 1 , 2 . 


el{X;ti,t2) = 
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Proof: Using by parts technique on the right hand side integral of fl3.12p . we get 


F{xi,t2 


\ogF{xi,t2)dxi = mf logF(ti,t2)- 


dF(xi,t2) 

dxi 


f*Xi 


F{u,t2) 


7o Fih,t2) - - ^ F{h,t2)\Jo F{x,,t2) 

Hence the hrst part follows from (13.12^ . The proof for second part is analogous. □ 


du 


On using Theorem 13.II and 13.21 we have the following result. The proof is omitted. 


Corollary 3.1 Forti,t 2 > 0, ifmf(ti,t 2 ) is increasing inti, i = 1,2 then e*{X-, 11 , 12 ) 
is also increasing in U. 


Let us set, for 0 ^ a < 6, 

Tl^\a,b-, t2) 



F{Xi,t2) 

F{ti,t2) 


dxi 


and 


4^\a,b-,ti) 



F{ti,X2) 
F{ti, ^ 2 ) 


dx2- 


Then an alternative expression to fl3.10l) and fl3.1ip is given hereafter. The proof is 
an immediate consequence of Fubini’s theorem. 


Theorem 3.3 Let X = {Xi,X 2 ) be a nonnegative bivariate random vector with 
e*{X-,ti,t 2 ) finite. Then 


U; 0 ) 


E 




Similarly, e 2 {X-,ti,t 2 ) 


E 


X2 ^(-^2,0;U) 1-^1 < ^5-^2 < ^2 


Let us discuss the effect of linear transformation on CDCPE. The proof is imme¬ 
diate from (13.1011 and (13.lip . 


Theorem 3.4 Let F and G be the bivariate distribution functions of two nonnegative 
random vectors X = {Xi,X 2 ) and Y = (Yi,U 2 ); respectively where Yi = ciXi + di 
with Ci > 0 and dj ^ 0 for i = 1,2. Then 


£iO^'^F,t2) = CiS* X; 


ti — di t2 — d2 


Cl C2 

which shows that CDCPE is a shift dependent measure. 


ti ^ di. 


Corollary 3.2 Let Y = (Id, 1 ^ 2 ) be a nonnegative bivariate random vector where 
Yi = CiXi + di with Ci > 0 and di ^ 0 for i = 1,2. Then, e^Y ; ti, ^ 2 ) is increasing in 
ti if and only z/£*(X; fi, ^ 2 ) is also so. 
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The effect of monotonic transformation on CDCPE is given in the next theorem. 

Theorem 3.5 Let X = {Xi,X 2 ) be a nonnegative bivariate random vector. Sup¬ 
pose ip{x) is strictly monotonic, continuous and differentiable function on [0, cxo) with 
</9(0) = 0 and derivative (p'{x). If ip is an increasing function with a ^ ip' ^ b, a,b > 0 
then 


be*{X;ip {ti),ip {t 2 )) ^e*{X^-,ti,t 2 ) ^ae*{X;ip {ti),ip (ts)), * = 1,2, 


where X^ = {p>{Xi), p>{X 2 )) is a bivariate random vector. If ip is decreasing with 
a ^ —ip' ^ b, then 

bei{X;ip~^{ti),ip~^{t2)) ^e*{X,^]ti,t2) ^ aei{X-,ip~^{ti),ip~^{t2)), 
where £j(X;ti,t 2 ), * = 1,2 is conditional DCRE as studied by Rajesh et al. (2014)- 
Proof: If ip is increasing then from fl3.10p and fl3.1ip . we have 

F{u,ip-\t2)) , F{u,p,-\t2)) 


^1’ ^2) — ~ 


and, el{X^-,ti,t 2 ) =- 


F {(p-^ih), (p-^{t2)) F {(p-^{ti), (p-^{t2)) 
F{cp-\R),v) , ^ Ficp-\R),v) 


log- 


ip'{u)du 
ip'{v)dv. 


F{(p ^iti),(p 1 (^ 2 )) F{(p ^{ti),(p 1 (^ 2 )) 

Hence the hrst part follows on using a ^ (p' ^ b. If 9 ? is decreasing we similarly obtain 
the second part of the proof. □ 


Theorem 3.6 For a strictly increasing function ipf) with (p(0) = 0, e*{X^;ti,t 2 ) is 
increasing in ti ife*{X]ti,t 2 ) is increasing in ti fori = 1,2 and p)' <1. 


Proof: Let p be strictly increasing. Then for i = 1, we get after some algebraic 
calculation 


d 


el{X^-ti,t2) 


sign 


dti ^{ti),p ^{t2)) 

rr~^ih) 

/ F{u,p~^{t 2 ))du / F{u, p~^{t 2 )) log F{u,p~^{t 2 ))p'{u)du 

Jo Jo 

- / F{u,p~^{t 2 ))p'{u)du / F{u,p~^{t 2 ))logF{u,p~^{t 2 ))d‘ 

Jo Jo 


^ 0 , 


where the last inequality follows on using that p' < 1. Thus, e\{X^p-,ti,t 2 ) is increasing 
in ti if £\{X-,p~^{ti),p~^{t 2 )), or equivalently, e*i{X-,ti,t 2 ) is increasing in R. The result 
follows analogously for i = 2. □ 


Now we define the following stochastic orders between two bivariate random vectors 
based on CDCPE. For more on stochastic order one may refer to Shaked and Shanthikumar 
(2007). 
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Definition 3.2 For two nonnegative bivariate random vectors X = (Xi,X 2 ) and Y = 
{Yi,Y 2 ), X is said to be greater than Y in CDCPE order (written as X ^cdcpe Y) if 
e*{X;ti,t 2 ) ^ e*{Y-,ti,t 2 ) for all 0, i = 1,2. 

Remark 3.2 It can be checked that the ordering defined above is reflexive, antisymmetric 
and transitive and thus partial ordering. 

Consider the following example to see that the ordering defined above is not equivalent 
to usual bivariate stochastic ordering. 

Example 3.3 Let X = {Xi,X 2 ) and Y = (Ei,!^) be two nonnegative bivariate random 
vector with distribution functions F (ti,t 2 ) and G {ti,t 2 ) respectively. Also let 

F{ti,t2) = -j_ -0 ^ ^ 4, 0 ^ ^ 4, 

with marginals Fx^iti) = ^, 0 ^ ^ 4 and Fx 2 (t 2 ) = 0 ^ t 2 ^ 4 and 

G{tl,t2) = ^ -j- 0 ^ ^ 1 , 0 ^ 12 ^ 1 , 

with marginals = tf, 0 ^ ^ 1 and 0 x 2 ( 12 ) = t^, 0 ^ t 2 ^ 1- From Figure\^ it 

is clear that for all 0 ^ ti ^ 1 and 0 ^ t 2 ^ 1 and for i = 1,2, £*(X\ti,t 2 ) ^ £*(Y; ti,t 2 }, 
giving X ^cdcpe Y. Again Fxi{0.4) — GYiiO.4) = —0.44 and Fxi(0.1) — Gyi (0.1) =0.015 
proving Xi '^st Yi. So, by Theorem 6.B.16 (c) of Shaked and Shanthikumar (2007) it can 
be concluded that X '^st Y. 



Graph of £l(X-,ti,t 2 ) - £\(Y-,ti,t 2 ) 


Graph of ^ 2 ) - e^iY] ti, ^ 2 ) 


Figure 2: Graphical representations of £*(X;ti, t 2 ) — £*(Y-, ti, ^ 2 ) (Example 13.3p 
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The following result can be obtained from Theorem 13.41 


Theorem 3.7 For two nonnegative bivariate random vectors X = (Xi,X 2 ) and X' = 
(X(, X 2 ), let Yi = CiXi + di and ¥( = CiX[ + di with Cj > 0 and di ^ 0 for i = 1,2. Then 
y ^CDCPE y if X ^cDCPE X' where Y' = {Yl,Yf). 

Now we have the more general result. 


Theorem 3.8 Let F,F',G and G' be the joint distribution functions of bivariate random 
vectors X, X', Y and Y', respectively. Also let for i = 1,2, Yi = OiXi + a and Yf = aiX[ + di 
with Oi > 0 and di ^ Ci > 0. Then, Y ^cdcpe Y' provided X ^cdcpe X' and either 
e*{X',ti,t 2 ) ore*{X';ti,t 2 ) is increasing in ti as well as t 2 . 


Proof: If X ^CDCPE X', then using Theorem 13.41 it can be written that 


i=*t\r \ ^ tl - Cl t2 - C2 

ei(Y-,ti,t2) ^ aiei[X-, -,- 

ai 02 

! ti — di t2 — d2 


> aiel[X'-, 


ai 


02 


= ri{Y'-,G,t2), 


where the last inequality follows on using the fact that e*{X'-,^, 12 ) is increasing in ti and 
t 2 . Similarly, the result follows if e*{X;ti,t 2 ) is increasing in ti and t 2 . □ 

In the recent past, the researchers have shown more interest in characterization of 
distributions. An important question regarding the CDCPE is whether it characterizes 
the underlying distribution function uniquely. In the following theorem we show that 
e*{X-,ti,t 2 ), i = 1,2 determines the distribution function uniquely. 


Theorem 3.9 Let X = (Xi, A 2 ) be a nonnegative bivariate random variable having ab¬ 
solutely continuous distribution function F with respect to the Lebesgue measure. Then 
CDCPE of X, defined in IS. 1 0\) and IS.111! , uniquely determines the distribution function 
provided they are finite. 


Proof: Let X and T be two bivariate random variables having joint distribution functions 
F and G, respectively. Also let, for all ti,t 2 ^ 0, 


e*{X-,ti,t 2 ) = e*{Y-,ti,t 2 ), i = l, 2. 


Differentiating e*{X-,ti,t 2 ) and e*{Y-,ti,t 2 ) with respect to U, i = 1,2 and, on using the 
relation cff {ti,t 2 )mf {ti,t 2 ) = 1 — -^mf {ti,t 2 ), we have from (I3.13P 


and 


■§f. {'mf{ti,t2)) 


-^e*{X-,ti,t2)mf{ti,t2) + e*{X-,ti,t2) -mf{ti,t2) 
e*{X-,ti,t2) - mf{ti,t2) 

{ti,t 2 ) + e*(T;ti,t 2 ) - rnj (ti,t 2 ) 
E*{Y-,ti,t 2 ) - mj{ti,t 2 ) 
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respectively. Suppose that 


e*{X;ti,t2) = e*{Y;tut2) = Oi (t) 


and 


for i 




(t) z + 9i{t)- z 
9i (t) - 2 ; 


1,2. Thus, we can write 


t = {tiM) , 


A 

dti 


{mf{ti,t2)) 




and 

— (mf (ti,t2)) = i’i (t,mf (ti,t2)) • 

In a recent paper, Thapliyal et al. (2013) proved that dynamic cumulative past entropy 
determines the distribution function uniquely, which in turn uniquely determines expected 
inactivity time. Thus e*{X-,ti,t 2 ) determines mf {ti,t 2 ), i = 1,2. Again, using the fact 
that vector valued expected inactivity time uniquely determines the bivariate distribution 
function, the proof is complete. □ 


Now we study the characterization result for uniform distribution. The proof follows 
easily as Xi and X 2 are independent. 

Theorem 3.10 Let X = (Ai, A 2 ) be a bivariate random variable with distribution function 
F. Then X is said to follow bivariate uniform distribution with distribution function 

F {ti,t2) = 0 ^ ^ 6, 0 i^t2 ^ d, 

ifandonlyife*{X;ti,t 2 ) = ^, i = l,2. 

In the following theorem we give another characterization result with dependent com¬ 
ponents. 

Theorem 3.11 Let X be a nonnegative bivariate random vector with e*{X;ti,t 2 ) finite 
and the components of bivariate FIT are mf {ti,t 2 ), i = 1,2 for all L ^ 0. Then, for 
0 < ti,t2 < 1, 0^ 0, 

£*(A;ti,t2) = i,j = l,2, i^j, (3.14) 

if and only if X is distributed as bivariate uniform with 

F{ti,t2) = 0 < ti,t2 < 1, 0. (3.15) 
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Proof: The if part is trivial on noting that if X follows the distribution (|3.15p then 




tj 

2 + 6 log tj ’ 


and £*{X]ti,t 2 ) 


ti{l + e\ogtj) . . , . 

(2 + 01ogt,)2’ 


i^j- 


To prove the converse let us assume that (I3.14h holds. Then, differentiating (I3.14h with 
respect to tj, and using (I3.13p . we get after some algebraic manipulation 


A 

dti 


mf{ti,t2) = 


1 


2 + 0 log tj 


hj = 1 , 2 , i^j 


which on integration gives 

mf{ti,t2) 


tj 

2 + 9 log tj 


+ Ci {tj ), 


where Ci{tj) is a constant of integration. Now, Ci{tj) = 0 as rnf{ti,t 2 ) —>■ 0 for ti —0, 
which in turn gives the bivariate EIT of (j3.15ll . Hence the result follows on using the fact 
that bivariate EIT determines the distribution uniquely. □ 


The following theorem gives a characterization of the bivariate power distribution. This 
result extends Theorem 6.2 of Di Crescenzo and Longobardi (2009) to bivariate setup. 

Theorem 3.12 Let X he a nonnegative bivariate random vector in the support (0, 6i) x 
( 0 , 62 ); h < 00 , i = l,2 with e*{X; ti, t 2 ) finite. Then 

£*{X-,ti,t2) = Ci{tj)mf{ti,t2), i,j = 1,2, ifi^j, (3.16) 


where Ci{tj) € ( 0 , 1 ) is a function independent of ti, characterizes the bivariate power dis¬ 
tribution 


F{ti,t2) 




6^0, 


(3.17) 


where Ci = cfibj)/ [1 - cfibj)]. 


Proof: The if part is straightforward, to prove the reverse implication, if (I3.16[) holds, then 
differentiating both the sides with respect to fi, and on using ()3.13p . we get after some 
algebraic manipulation 

d 

—(ti,t2 ) = [1 - Ci(tj)] , i,j = 1,2, i^j 
which on integration gives 


mf{ti,t2) = [1 - Ci{tj)] ti + ki{tj), 

where ki{tj) is a constant of integration. Now, ki{tj) = 0 as rnf{ti,t 2 ) —>■ 0 for ti —>■ 0. The 
rest of the proof follows from Theorem 2.1 of Nair and Asha (2008). 
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3.2 Conditional dynamic CPE for Xi given Xj = tj 


The determination of the joint distribution function of X = {Xi,X 2 ), when conditional 
distributions of {Xi\X 2 = ^ 2 ) and {X 2 \Xi = ti) are known, has been an important problem 
dealt with by many researchers in the past. This approach of identifying a bivariate density 
using the conditionals is called the conditional specification of the joint distribution (see 
Arnold et ah, 1999). These conditional models are often useful in many two component 
reliability systems, when the operational status of one component is known. Let the dis- 
tribution function of Yi = {Xi\Xi < ti,Xj = tj), i,j = 1,2, i ^ j be F* {ti\tj). Then, for 
an absolutely continuous nonnegative bivariate random vector X, the conditional dynamic 
CPE of Yi is defined as 


-f*{X;F,t2) = - 


f 


Fi {xi\tj) T) 


F'iitiltj 


log 


F'iitiltj 


-dxi, 


Xi < ti 


(3.18) 


i,j = 1,2, i jtz j. In particular, if Xi and X 2 are independent, then (j3.18p reduces to 
marginal dynamic CPE of Xi, i = 1,2 as given in (11.51) . Eollowing Roy (2002) the bi- 

—X 

variate reversed hazard rate of X = {Xi,X 2 ) is also defined by a vector, cf) {ti\tj) = 
(ii|i 2 ),?f (i 2 |ii)), where {ti\tj) = J-logT)* (ti|tj), i,j = 1,2, i / j. For i = 1, 

(^ 11 ^ 2 )^^! is the probability of failure of the first component in the interval {ti — Ati, ti] 
given that it has failed before ti and the failure time of the second is t 2 - Another defini¬ 
tion of bivariate EIT of X = (Xi, A 2 ) is given by Kayid (2006) as a vector, {ti\tj) = 
{m^{ti\t 2 ),rn^{t 2 \ti)), where rrif {ti\tj) = E{ti-Xi\Xi <ti,Xj =tj), i,j = 1 , 2 , i / j. 
For i = 1, 

mf{ti\t2) = FI {xi\t2) dxi, 

which measures the expected waiting time of Xi given that Xi < ti and X 2 = t 2 - Unlike 
4>^{ti,t2) and {ti,t 2 ), fn^ {ti\tj) determines the distribution uniquely. But, (j)^ {ti\tj) 
does not provide F{ti,t 2 ) uniquely. 

Differentiating (|3.18|) with respect to ti and simplifying, we get 


A 

dti 


7 *(X;ti,t2) = A {ti\tj) [mf {ti\tj) - 7-(X;ti,t2)] , i,j = l,2,z ^ j. 


Now we have the following theorem. 


Theorem 3.13 Forti,t 2 > 0, 7 *(X;ti,t 2 ) is increasing inti, if and only if 

li{X]ti,t2) ^mf {ti\tj), i,j = l,2,i ^ j. 


The following example gives an application of the above theorem. 


Example 3.4 If the density function of a continuous bivariate random vector X = (Xi,X 2 ), 
/ (a^i) ^ 2 ) is given by 


f{xi,X2) 


J{xi + 4 x 2 ), 0 ^ xi ^ 2, 0 ^ X2 ^ 1 

0 , elsewhere. 
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then it can he checked that for i,j = 1,2 and i ^ j, j*{X;ti,t 2 ) is increasing in U for all 
tj ^ 0. Again, FigurelE shows that for 0 ^ ti ^ 2 and 0 ^ t 2 ^ 1 and i = 1,2, i ^ j, 
mf-— 7 *(X;ti,t 2 ) are always positive, satisfying Theorem \3.13[ 




Graph of mf (ti|t 2 ) - 

Figure 3: Graph of mf 


Graph of (^ 2 !^) - 72 (-^;GG 2 ) 
7 *(X;fi,f 2 ) fExample I3.4p 


The next theorem, analogous to Theorem I3.21 establishes relation between dynamic 
CPE and EIT of Xi given Xj = tj, i,j = 1, 2, i 7 ^ j. The proof is omitted. 

Theorem 3.14 For ti,t 2 > 0, let ^*{X-,ti,t 2 ) andnif{ti\tj) are CDCPE and the com¬ 
ponents of bivariate EIT as defined above. If /* {ti\tj) denotes the density function of 
{Xi\Xi < ti,Xj =tj), then 


r 

li{X-,ti,t 2 ) = / nif {xi\tj) f* {xi\tj)dxi, i,j 
Jo 


= 1 , 2 , i^j. 


Let f {ti,t 2 ) be the joint pdf of the bivariate random variable X, then (I3.18P can 
alternatively be written as 


lliX'MM) = {ti\tj)\og f f{z,tj)dz- [ 

Jo Jo 

hi = 1 , 2 , i j. So, taking 




log / f{z,tj)dzds, 


Ff {ti\tj) Jo 


Tf\a,b) = - log 


[ f{z,tj) 

Jo 


dz ds, i,j = 1,2, i / j. 


and proceeding in the similar way as of previous subsection, it can be shown that 


'y*{X;ti,t 2 ) = mf {ti\tj) log I f{z,tj)dz+E TY\Xi,ti)\Xi < ti,Xj = tj , i, j = 1,2, i ^ j. 


( 2 ), 


18 





















Again, for 0 ^ a ^ 6, defining 


^ (a, b) = - I log 


fo f tj) dz 


ds, i,j = 1,2, j, 


° !q f {^^'tj)dz 

as before, we have the following theorem, which is analogous to Theorem 13.31 


Theorem 3.15 LetX = {Xi,X 2 ) be anonnegative bivariate random vector with j*{X]ti,t 2 ) < 
oo, i = 1,2. Then 


ri{x-,h,t2) = E 


T\^\Xi,ti)\Xi < ti,Xj = tj 


i,j = 1,2, i / j. 


The effect of linear transformation on ^*{X;t\,t 2 ) is given in the following theorem. 


Theorem 3.16 Let X = (Xi, A 2 ) and Y = (Yi,!^) be two nonnegative bivariate random 
variables having distribution functions F and G respectively, where Yi = CiXi + di, i = 1,2. 
Then 


r^{Y■M,t2) = cir^[x 


ti — di t2 — d2 


Cl 


C2 


ti ^ di. 


Corollary 3.3 Let Y = {Yi,Y 2 ) he a nonnegative bivariate random vector where Yi = 
CiXi + di with Ci > 0 and di ^ Q for i = 1,2. Then, 7*(Y; ti, ^ 2 ) is increasing in ti if and 
only if {X; ti, t 2 ) is also so. 


The immediate consequence of Theorem 13.161 are the following two theorems. The proof 
being analogous to Theorems 13.71 and 13.81 is omitted. 

Theorem 3.17 LetX = {Xi,X2) and X' = (X[,X2) be two nonnegative bivariate random 
variables. Also let Y = (Yi,Y2) and Y' = {Yl,Yf) he two nonnegative bivariate random 
vectors such that Yi = CiXi + di and Yf = CiX[ + di with Ci > 0, di ^ 0 for i = 1,2. 
If ^*{X;ti,t2) ^ j*{X';ti,t2) then j*{Y;ti,t2) ^ ^*{Y'-,ti,t2), for all ti,t2 ^ 0 , where 
F,F',G,G' are the distribution functions of X, X',Y andY', respectively. 


Theorem 3.18 Let F, F',G,G' are the distribution functions of nonnegative bivariate ran¬ 
dom variables X, X', Y and Y', respectively. Also let for i = 1,2, Yi = OiXi + Ci and Yf = 
aiXl -\-di with a, > 0 and di ^ Ci > 0. Then, for all ti,t2 ^ 0, j*{Y-,ti,t2) ^ j*{Y'-,ti,t2) 
provided ^*{X;ti,t2) ^ 7*(X';ti,t2) and either {X ■, ti, t2) or ^*{X'-,ti,t2) is increasing 
in ti as well as ^2- 
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